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Abstract 



In a previous paper we analyzed fermion masses (focusing on neu- 
trino masses and mixing angles) in an SO(10) SUSY GUT with 
U(2)xU(l) ?l family symmetry. The model is "natural" contain- 
ing all operators in the Lagrangian consistent with the states and 
their charges. With minimal family symmetry breaking vevs the 
model is also predictive giving a unique solution to atmospheric 
(with maximal — » v T mixing) and solar (with SMA MSW 
v e — > v s mixing) neutrino oscillations. In this paper we analyze 
the case of general family breaking vevs. We now find several new 
solutions for three, four and five neutrinos. For three neutrinos 
we now obtain SMA MSW, LMA MSW or vacuum oscillation 
solutions for solar neutrinos. In all three cases the atmospheric 



data is described by maximal z/ M — > v T mixing. In the four and 
five neutrino cases, in addition to fitting atmospheric and solar 
data as before, we are now able to fit LSND data. All this is 
obtained with the additional parameters coming from the fam- 
ily symmetry breaking vevs; providing only minor changes in the 
charged fermion fits. 



+ On leave of absence from Faculty of Mathematics and Physics, Comenius 
Univ., Bratislava, Slovakia 
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1 Introduction 



Neutrino oscillations provide a window onto new physics beyond the standard 
model Q and several experiments now provide evidence for neutrino oscilla- 
tions. This includes data on solar neutrinos 0, atmospheric neutrinos 
and the accelerator-based experiment, liquid scintillator neutrino detector 
[LSND] ||. These positive indications are constrained by null experiments 
such as Chooz |4j] and Karmen |J. The data strongly suggests that neutrinos 
have small masses and non- vanishing mixing angles In the near future, 
many more experiments will test the hypothesis of neutrino masses fj, - 
[|T3j]. Thus there is great excitement and anticipation in this field. 

In a recent paper I |14| we analyzed an SO w super symmetric [SUSY] 
grand unified theory [GUT] with family symmetry Ui x [7f . The theory 
was "natural," i.e. the Lagrangian was the most general consistent with the 
states and symmetries. In addition, with minimal family symmetry break- 
ing vacuum expectation values [vevs], the number of arbitrary parameters 
in the effective low energy theory, below the GUT scale, was less than the 
number of observables. Hence the theory was "predictive" and testable. We 
analyzed the predictions for charged fermion masses and mixing angles us- 
ing a global x 2 analysis [T5|, TA\ finding excellent agreement with the data. 
In the neutrino sector we obtained a unique solution to both atmospheric 
and solar |l]] neutrino oscillation data. This solution has three active and one 
sterile neutrino. It has maximal — > v T oscillations fitting atmospheric data 
and small mixing angle [SMA] Mikheyev, Smirnov, Wolfenstein [MSW] [[L6[ 
v e v s (where s denotes sterile) oscillations for solar data, without fine- 
tuning. We were however unable to simultaneously fit LSND ||, even with 
four neutrinos. In addition, we were unable to find a three neutrino solution 
to both atmospheric and solar neutrino data. It is imperative to understand 
if these results are robust. In particular, without a theory of family sym- 
metry breaking we may consider more general family symmetry breaking 
vevs.Q In this paper we allow for the most general family symmetry break- 
ing vevs; introducing two new complex parameters /t(i,2)- There are now 

1 In the Standard Model, the three active neutrino species (members of electroweak 
doublets) are massless. As a consequence individual lepton number is conserved and 
neutrinos cannot oscillate. 

2 In fact, we noted in [|l4| that it is possible to obtain three neutrino solutions to 
atmospheric and solar data if we allow for non-minimal family symmetry breaking vevs. 
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more parameters for charged fermion masses and mixing angles than there 
are observables. The new parameters have minor consequences for charged 
fermions (fits to m e , m^, and V us , which are all known to excellent accuracy, 
require them to remain small), but significant consequences for neutrinos. In 
fact, with the additional parameters we are now able to obtain three possible 
three-neutrino solutions to atmospheric || and solar [|IJ neutrino data. With 
one or two sterile neutrinos we can also obtain solutions to atmospheric [§], 
solar and LSND @ data. 

In section 2, we discuss the model and family symmetry breaking. The 
model is an SO(10) [SUSY GUT] xU(2)xU(l) n [family symmetry] model. It 
is a small variation of the theory introduced by Barbieri et al. [BHRR] |17| 
where the non-abelian family symmetry was introduced to provide a natural 
solution to flavor violation in SUSY theories [[R| [19], |2(J. In section 3, we 



present the general framework for neutrino masses and mixing angles. In 
section 4, we describe the three neutrino solutions and in sections 5 and 6 
we present the four and five neutrino solutions, respectively. Our conclusions 
are in section 7. 



2 An SO(10)xU(2)xU(l)x model 

The three families of fermions are contained in 16 a , a = 1,2; and I63 where 
a is a U(2) flavor index. [Note U(2) = SU(2) x U(l)' where the U(l)' charge 
is +1 (—1) for each upper (lower) SU(2) index.] At tree level, the third 
family of fermions couples to a 10 of Higgs with coupling A I63 10 I63 in the 
superspace potential. The Higgs and I63 have zero charge under the first 
two U(l)s, while 16 a has charge —1 and thus does not couple to the Higgs 
at tree level, f] 

3 There are in fact four additional U(l)s implicit in the superspace potential (eqn. ||). 
These are a Peccei-Quinn symmetry in which all 16s have charge +1, all 16s have charge 
— 1, and 10 has charge —2; a flavon symmetry in which (<f> a , S ab , A ab ) and M have charge 
+1 and \b has charge —1; a symmetry in which M\ M" have charge +1 and x, X a have 
charge -1 and and an R symmetry in which all chiral superfields have charge +1. The 
family symmetries of the theory may be realized as either global or local symmetries. For 
the purposes of this paper, it is not necessary to specify which one. However, if it is 
realized locally, as might be expected from string theory, then not all of the U(l)s are 
anomaly free. We would then need to specify the complete set of anomaly free U(l)s. 
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Three superfields (0°, S ab = S ba , A ab = — A ba ) are introduced to spon- 
taneously break U(2)xU(l) and to generate Yukawa terms giving mass to 
the first and second generations. The fields (0°, S ab , A ab ) are SO(10) sin- 
glets with U(l) charges {0, 1, 2}, respectively. The most general vacuum 
expectation values are given by 

(0 2 > ± o, 

(s 22 )^o, (S") = Kl (S"), <s 12 > = ms 22 >, (i) 

(A 12 ) + 0, 

where the constants K\, k<i are arbitrary. The vevs ((0 2 ) ~ (S 22 ) ~ tM^/ (45)) 
break U(2)xU(l) to U(l) and (A 12 ~ e'M ) completely. In this model, sec- 
ond generation masses are of order e, while first generation masses are of order 
e /2 /e. In paper I HU we analyzed this theory with minimal family breaking 
vevs (ki = K2 — 0). In this paper we show the effects of non- vanishing «(i.2)- 
The superspace potential for the charged fermion sector of this theory, 



including the heavy Froggatt-Nielsen states [21], is given by 



WD I63IOI63 + 16 a 10x a (2) 

+ Xa (M X a + <P a X + S ab X b + A ab 16 6 ) 
+ x a (M' Xa + 45 16 a ) 
+ x (M" X + 45 16 3 ) 

where 

M = M (1 + a X + PqY). (3) 

X, Y are SO(10) breaking vevs in the adjoint representation with X corre- 
sponding to the U(l) in SO(10) which preserves SU(5), Y is standard weak 
hypercharge and ao, Po are arbitrary parameters. The field 45 is assumed to 
obtain a vev in the B—L direction. Note, this theory differs from [BHRR] P"Tf 
only in that the fields <p a and S ab are now SO(10) singlets (rather than SO(10) 
adjoints) and the SO(10) adjoint quantum numbers of these fields, necessary 
for acceptable masses and mixing angles, has been made explicit in the field 
45 with U(l) charge l|] This theory thus requires much fewer SO (10) ad- 
joints. Moreover our neutrino mass solution depends heavily on this change. 

4 This change (see BHRR |i~7|]) is the reason for the additional U(l)s. 
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The effective mass parameters Mo, M', M" are SO(IO) invariants.^ The 
scales are assumed to satisfy M ~ M' ~ M" > (0 2 ) ~ (S 22 ) > (A 12 ) 
where M may be of order the GUT scale. In the effective theory below 
M , the Froggatt-Nielsen states {x, X, X° '■> Xa, Xa, X a } ma y De integrated 
out, resulting in the effective Yukawa matrices for up quarks, down quarks, 
charged leptons and the Dirac neutrino Yukawa matrix given by (see fig. 1) [] 
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(5) 



jy = (Baryon^ — Lepton#) 
for f — {Q,u,d,L,e,u}. 



(6) 



5 The effective mass parameters represent vevs of SO(10) singlet chiral superfields. 

6 Note, we use the notation of BHRR fl7| . The parameter p vanishes in the limit 
(h = 0(see equations ^). This is a consequence of the B-L vev in the 2-2 entry or 
the anti-symmetry of the coupling to A ab in the 1-2 clement which is in conflict with the 
SU(5) invariance of M in this limit which only allows for symmetric u — u couplings. 
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Figure 1: Diagrams generating the Yukawa matrices 

In our notation, fermion doublets are on the left and singlets are on the 
right. Note, we have assumed that the Higgs doublets of the minimal super- 
symmetric standard model[MSSM] are contained in a single 10 dimensional 
SO(10) multiplet. Hence all the fits have large values of tan/?. [] 

2.1 Results for Charged Fermion Masses and Mixing 
Angles 

We have performed a global x 2 analysis, incorporating two (one) loop renor- 
malization group [RG] running of dimensionless (dimensionful) parameters 
from Mq to Mz in the MSSM, one loop radiative threshold corrections at 
M z , and 3 loop QCD (1 loop QED) RG running below M z . Electroweak 
symmetry breaking is obtained self-consistently from the effective potential 
at one loop, with all one loop threshold corrections included. This analysis 

7 Note, we could obtain small values of tan /3 in SO(10) at the cost of one new parameter. 
If the 10 which couples to fermions mixes with other states then the Higgs field coupling 
to up and down quarks may have different effective couplings to matter, i.e. such that 
A 10 D A H u + £ Hd- We could then consider two limits — case (1) A = £ (no Higgs 
mixing) with large tan/3, and case (2) A £ or small tan/3. In paper I, we also considered 
case (2) and found no significant improvements in the fit. 
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Table 1: Charged fermion masses and mixing angles 

Initial parameters: K\ — K2 — 

(l/a G , Mq, e 3 ) = (24.52, 3.03 • 10 16 GeV, -4.06%), 

(A,r, a, e, p, e') = (0.79, 12.4, 0.84, 0.011, 0.043, 0.0031), 

($ CT , $ e , $ p ) = (0.73, -1.21, 3.72)rad, 

(m , Mi /2 , A), /i(M z )) = (1000, 300, -1431, 110) GeV, 

((m^/mo) 2 , (m Hu /m ) 2 , tan/5) = (2.23, 1.66, 53.7) 



Observable 


Data(cr) 


Theory 




(masses 


in GeV) 


M z 


91.187 (0.091) 


91.17 


M w 


80.388 (0.080) 


80.39 


G> • 10 5 

a EM 

a s (M z ) 
■ 10 3 


1.1664 (0.0012) 
137.04 (0.14) 
0.1190 (0.003) 
-0.20 (1.1) 


1.166 
137.0 
0.1174 
+0.314 


M t 


173.8 (5.0) 


174.9 


m b (M b ) 
M b -M c 


4.260 (0.11) 
3.400 (0.2) 


4.331 
3.426 


m s 

m d /m s 

Q- 2 

M T 
M, 

M e ■ 10 3 


0.180 (0.050) 
0.050 (0.015) 
0.00203 (0.00020) 

1.777 (0.0018) 
0.10566 (0.00011) 
0.5110 (0.00051) 


0.147 

0.0589 

0.00201 

1.777 

.1057 

0.5110 


v us 


0.2205 (0.0026) 


0.2205 


v cb 


0.03920 (0.0030) 


0.0403 


v ub /v cb 


0.0800 (0.02) 


0.0691 


B K 


0.860 (0.08) 


0.870 


B{b^s-i)-W 


3.000 (0.47) 


2.992 


TOTAL x 2 


2.26 
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is performed using the code of Blazek et.al. ]T3|. Q In this paper, we just 



present the results for one set of soft SUSY breaking parameters mo, M\m 
with all other parameters varied to obtain the best fit solution. In the first 
two columns of table [I] we give the 20 observables which enter the \ 2 func- 
tion, their experimental values and the uncertainty a (in parentheses). In 
most cases a is determined by the 1 standard deviation experimental uncer- 
tainty, however in some cases the theoretical uncertainty (~ 0.1%) inherent in 
our renormalization group running and one loop threshold corrections dom- 



inates. Lastly, in contrast to paper I we include a 1999 updated value p2 
of Pnew, the measure of SU(2) violation beyond the standard model. This 
change substantially improves our global charged fermion fits. 

There are 8 real Yukawa parameters and 5 complex phases. We take the 
complex phases to be <3> p , $ e , § a , $ Kl and $ K2 . With 13 fermion mass 
observables (charged fermion masses and mixing angles replacing €k as 
a "measure of CP violation" ]) we have enough parameters to fit the data. 
In table [I] we also show the fits obtained with K\ = k 2 = as a benchmark 
for the cases with non-zero Ki^ which follow. From table [I] it is clear that 
this theory fits the low energy data quite well. Q 

Finally, the squark, slepton, Higgs and gaugino spectrum of our theory is 
consistent with all available data. The lightest chargino and neutralino are 
higgsino-like with the masses close to their respective experimental limits. 
As an example of the additional predictions of this theory consider the CP 
violating mixing angles which may soon be observed at B factories. For the 
selected fit with K\ = k 2 = we find 

(sin 2a, sin 2(3, sin 7) = (0.74,0.54,0.99) (7) 



8 We assume universal scalar mass mo for squarks and sleptons at Mq. We have not 
considered the flavor violating effects of U(2) breaking scalar masses in this paper. 

9 The Jarlskog parameter J = Im(V u dV* b V c bV* d ) is a measure of CP violation. We 
test J by a comparison to the experimental value extracted from the well-known K° — K° 
mixing observable £k = (2.26±0.02) x 10~ 3 . The largest uncertainty in such a comparison, 
however, comes in the value of the QCD bag constant Bk ■ We thus exchange the Jarlskog 
parameter J for Bk in the list of low-energy data we are fitting. Our theoretical value 
of Bk is defined as that value needed to agree with ck for a set of fermion masses and 
mixing angles derived from the GUT-scale. We test this theoretical value against the 
"experimental" value of Bk- This value, together with its error estimate, is obtained from 



recent lattice calculations 23 . 

10 Note, the strange quark mass m s (lGeV) ~ 150 MeV is small, consistent with recent 
lattice results. 
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or equivalently the Wolfenstein parameters 



(p, 77)= (-0.04,0.31) . (8) 

As an aside, we have also computed the SUSY contribution to the muon 
anomalous magnetic moment. Our prediction for the selected SUSY point [^] 
gives values for a^ USY ~ 40 • 10 -10 , in good agreement with the latest pre- 



liminary data from the ongoing BNL experiment |25 

In tables | - | we give results for non-zero K\, K2- These results have 
been obtained with a slightly different procedure than previously. We have 
followed a multi step iterative procedure for finding "good" fits to both 
charged fermion and neutrino data. This is in lieu of combining the neu- 
trino and charged fermion sectors into a single y 2 function and minimizing 
the total x 2 with respect to variations of all the parameters. Let us now 
describe this procedure in more detail. 

In each case we select a pair of non-zero values for k\ and k,2 and keep 
these two parameters fixed while we repeat the charged fermion analysis. 
If we obtain a good fit, we use these as initial values for the analysis of 
the neutrino sector (discussed in the next section). Then in the neutrino 
analysis we only vary those parameters not already included in the charged 
fermion analysis. If the resulting neutrino fit is not acceptable, we make a 
step in the (k\, Kz) parameter space and start again with the charged fermion 
analysis. We also found that we can improve the neutrino fit for fixed K\ and 
K2 if we return to the charged fermion analysis and carefully move one or 
more parameters entering the Yukawa matrices slightly away from their best 
fit value (watching so as not to incur large changes in the charged fermion 
contributions to y 2 )- Thus our tables ^|- |5| do not show the absolute "best" 
fits for fixed K\ and K2. Following this procedure we focus independently on 
different neutrino solutions as indicated in the table captions. Thus although 
the data in tables |2]-|^ do not seem much different, they do however represent 
significant changes in the neutrino sector; discussed in the next section. 

Before we conclude this section, let us consider one test (in the charged 
fermion sector) which may be able to distinguish among these different neu- 
trino cases. The unitarity angles (sin 2a, sin 2f3, sin 7) or equivalently the 

11 Although this result does depend on the particular point in SUSY parameter space 
we have selected, it is independent of the particular neutrino solution. In addition, we 
have assumed universal masses for squarks and sleptons at the GUT scale. Non-universal 
slepton masses can affect our result. 
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Table 2: Charged fermion masses and mixing angles: 3 neutrino 
SMA MSW 

Initial parameters: K\ — k%, \k 2 \ — 0.028 

(l/a G , Mq, e 3 ) = (24.52, 3.05 • 10 16 GeV,-4.07%), 

(A,r, a, e, p, e') = (0.79, 12.3, -0.96, 0.010, 0.042, 0.0031), 

($ CT , $ e , $ p , $ K1 , $ K2 ) = (3.84, 0.0032, 5.02, -1.70, -0.85)rad, 

(mo, Mi/a, A), /i(M z )) = (1000, 300, -1438, 110) GeV, 

((m Hd /™o) 2 , (m Hu /m ) 2 , tan/5) = (2.22, 1.66, 53.7) 



vy uoci vcxuic 


i^J CI L d I C 1 


/ h pnvn 

J. iLy^WI U 






in GeV) 

111 v_ J V_ V 1 


My 


91 187 fO 091) 


91.18 


M w 


80.388 (0.080) 


80.40 


G> • 10 5 

a EM 

a s (M z ) 
o ■ 10 3 


1.1664 (0.0012) 
137.04 (0.14) 
0.1190 (0.003) 
-0.20 (1.1) 


1.166 
137.0 
0.1174 
+0.322 


M t 


173.8 (5.0) 


175.0 


m 6 (M ft ) 
M b -M c 


4.260 (0.11) 
3.400 (0.2) 


4.326 
3.432 


m s 

m d /m s 

Q- 2 

M T 


0.180 (0.050) 
0.050 (0.015) 
0.00203 (0.00020) 
1.777 (0.0018) 


0.146 
0.0585 
0.00201 
1.776 


M, 


0.10566 (0.00011) 


.1057 


M e ■ 10 3 

V us 

v ub /v cb 


0.5110 (0.00051) 
0.2205 (0.0026) 
0.03920 (0.0030) 
0.0800 (0.02) 


0.5110 
0.2206 
0.0402 
0.0702 


B K 


0.860 (0.08) 


0.8691 


B(b^s-i)-W 


3.000 (0.47) 


2.958 


TOTAL x 2 


2.48 
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Table 3: Charged fermion masses and mixing angles: 3 neutrino 
LMA MSW 

Initial parameters: |«i| = 0.055, \k 2 \ — 0.31 

(l/a G , Mq, e 3 ) = (24.52, 3.05 • 10 16 GeV, -4.08%), 

(A,r, a, e, p, e') = (0.79, 14.3, -1.13, 0.009, 0.045, 0.0028), 

($ CT , $ e , $ p , $ K1 , $ K2 ) = (3.82, -0.69, 4.83, 4.07, -1.14)rad, 

(m , Mi /2 , A , ii(M z )) = (1000, 300, -1444, 110) GeV, 

{{m Hd /m )\ (m H Jm ) 2 ,t&n(3) = (2.22, 1.66, 53.7) 



vy uoci vciuic 


i^J CI L d I C J 


/ h pnvn 

J. iLy^WI U 






in GeV) 

111 v_ J V_ V 1 


M 7 


91 187 fO 091) 


91.18 


M w 


80.388 (0.080) 


80.40 


G> • 10 5 

a s (M z ) 
o ■ 10 3 


1.1664 (0.0012) 
137.04 (0.14) 
0.1190 (0.003) 
-0.20 (1.1) 


1.166 
137.0 
0.1174 
+0.322 


M t 


173.8 (5.0) 


174.9 


m 6 (M ft ) 
M b -M c 


4.260 (0.11) 
3.400 (0.2) 


4.323 
3.433 


m s 

m d /m s 

Q- 2 

M T 


0.180 (0.050) 
0.050 (0.015) 
0.00203 (0.00020) 
1.777 (0.0018) 


0.138 
0.0664 
0.00202 
1.776 


M, 


0.10566 (0.00011) 


.1057 


M e ■ 10 3 

V us 

v ub /v cb 


0.5110 (0.00051) 
0.2205 (0.0026) 
0.03920 (0.0030) 
0.0800 (0.02) 


0.5110 
0.2204 
0.0409 
0.0782 


B K 


0.860 (0.08) 


0.8682 


B(b^s-i)-W 


3.000 (0.47) 


2.999 


TOTAL x 2 


3.99 
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Table 4: Charged fermion masses and mixing angles: 3 neutrino 
Vacuum 

Initial parameters: |«i| = 0.004, | /« 2 1 = 0.025 

(l/a G , Mq, e 3 ) = (24.52, 3.05 • 10 16 GeV, -4.16%), 

(A,r, a, e, p, e') = (0.80, 15.6, -0.35, 0.013, 0.041, 0.0035), 

($ CT , $ e , $ p , $ K1 , $ K2 ) = (3.00, -0.65, 4.41, 3.74, -0.052)rad, 

(m , Mi /2 , A , fx(M z )) = (1000, 300, -1433, 110) GeV, 

{{m Hd /m )\ (m H /m ) 2 , tan/3) = (2.22, 1.66, 53.7) 



vy uoci vciuic 


i^J CI L d I C J 


/ h pnvn 

J. iLy^WI U 






in GeV) 

111 v_ J V_ V 1 


M 7 


91 187 fO 091) 


91.18 


M w 


80.388 (0.080) 


80.40 


G> • 10 5 

a s (M z ) 
o ■ 10 3 


1.1664 (0.0012) 
137.04 (0.14) 
0.1190 (0.003) 
-0.20 (1.1) 


1.166 
137.0 
0.1171 
+0.322 


M t 


173.8 (5.0) 


175.0 


m 6 (M ft ) 
M b -M c 


4.260 (0.11) 
3.400 (0.2) 


4.324 
3.405 


m s 

m d /m s 

Q- 2 

M T 


0.180 (0.050) 
0.050 (0.015) 
0.00203 (0.00020) 
1.777 (0.0018) 


0.170 
0.0548 
0.00202 
1.776 


M, 


0.10566 (0.00011) 


.1057 


M e ■ 10 3 

v us 
v cb 

v ub /v cb 


0.5110 (0.00051) 
0.2205 (0.0026) 
0.03920 (0.0030) 
0.0800 (0.02) 


0.5110 
0.2205 
0.0392 
0.0758 


B K 


0.860 (0.08) 


0.8604 


B(b^s-i)-W 


3.000 (0.47) 


2.938 


TOTAL x 2 


1.47 
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Table 5: Charged fermion masses and mixing angles: 4 neutrino 
SMA MSW + LSND 

Initial parameters: \k,i\ — 0.0001, \k 2 \ — 0.002 
(l/a G , Mq, e 3 ) = (24.50, 3.07- 10 16 GeV,-4.14%), 
(A,r, a, e, p, e') = (0.75, 12.4, -0.76, 0.011, 0.044, 0.0032), 
($ CT , $ e , $ p , $ K1 , $ K2 ) = (3.87, -0.95, 3.97, 4.81, 1.13)rad, 
(m , Mi/a, 4,, /i(M z )) = (1000, 300, -1459, 110) GeV, 
((m^/™o) 2 , ( mH „/m ) 2 , tan/5) = (2.19, 1.65, 53.0) 



( ) neiP'FVP nip 


i^J CI L d I C J 


/ h pnvn 

J. iLy^WI U 






in GeV) 

111 v_ < V_ V 1 


M 7 


91 187 fO 091) 


91.18 


M w 


80.388 (0.080) 


80.40 


G> • 10 5 

a s (M z ) 
o ■ 10 3 


1.1664 (0.0012) 
137.04 (0.14) 
0.1190 (0.003) 
-0.20 (1.1) 


1.166 
137.0 
0.1173 
+0.318 


M t 


173.8 (5.0) 


173.5 


m 6 (M ft ) 
M b -M c 


4.260 (0.11) 
3.400 (0.2) 


4.341 
3.422 


m s 

m d /m s 

Q- 2 

M T 


0.180 (0.050) 
0.050 (0.015) 
0.00203 (0.00020) 
1.777 (0.0018) 


0.148 
0.0591 
0.00201 
1.776 


M, 


0.10566 (0.00011) 


.1057 


M e ■ 10 3 

V us 

v ub /v cb 


0.5110 (0.00051) 
0.2205 (0.0026) 
0.03920 (0.0030) 
0.0800 (0.02) 


0.5110 
0.2205 
0.0402 
0.0699 


B K 


0.860 (0.08) 


0.8696 


B(b^s-i)-W 


3.000 (0.47) 


3.007 


TOTAL x 2 


2.94 
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Table 6: Charged fermion masses and mixing angles: 5 neutrino 
SMA MSW + LSND 

Initial parameters: |«i| = | ac 2 1 2 , \^2\ = 0.032 

(l/a G , M G , e 3 ) = (24.52, 3.06 • 10 16 GeV,-4.09%), 

(A,r, a, e, p, e') = (0.79, 12.2, -0.94, 0.011, 0.042, 0.0031), 

($ CT , $ e , $ p , $ K1 , $ K2 ) = (3.84, 0.07, 5.03, -2.49, -1.19)rad, 

(m , Mi/a, 4,, /i(M z )) = (1000, 300, -1438, 110) GeV, 

((m Hd /m ) 2 , (m H Jm ) 2 ,tanp) = (2.22, 1.66, 53.7) 



vy uoci vciuic 


i^J CI L d I C J 


/ h pnvn 

J. liOf/ U 






in GeV) 

111 VJv V 1 


My 


91 187 fO 091) 


91.17 


M w 


80.388 (0.080) 


80.40 


G> • 10 5 

a EM 

a s (M z ) 
o ■ 10 3 


1.1664 (0.0012) 
137.04 (0.14) 
0.1190 (0.003) 
-0.20 (1.1) 


1.166 
137.0 
0.1174 
+0.322 


M t 


173.8 (5.0) 


175.0 


m b (M b ) 
M b -M c 


4.260 (0.11) 
3.400 (0.2) 


4.328 
3.426 


m s 

m d /m s 

Q- 2 

M T 


0.180 (0.050) 
0.050 (0.015) 
0.00203 (0.00020) 
1.777 (0.0018) 


0.148 
0.0588 
0.00201 
1.777 


M, 


0.10566 (0.00011) 


.1057 


M e ■ 10 3 

V us 

v ub /v cb 


0.5110 (0.00051) 
0.2205 (0.0026) 
0.03920 (0.0030) 
0.0800 (0.02) 


0.5110 
0.2205 
0.0401 
0.0701 


B K 


0.860 (0.08) 


0.8686 


B(b^s-i)-W 


3.000 (0.47) 


2.983 


TOTAL x 2 


2.12 
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Table 7: Unitarity triangle angles and Wolfenstein parameters for 
the different neutrino fits with non zero K\, k 2 - 



Neutrino fit 


Values of 


K X , 


^2 




(sin 2a, sin 2/5, sin 7) 


(p, v) 


3i> SMAMSW 


Kx 




K 2 \ 


= 0.028 


(0.92, 0.50, 0.95) 


(-0.10, 0.30) 


3u LMAMSW 




Ki 


= 0.055, 




= 0.31 


(0.94, 0.39, 0.73) 


(-0.24, 0.26) 


3u Vacuum 




K! 


= 0.004, 


^2 


= 0.025 


(0.86, 0.56, 0.97) 


(-0.08, 0.33) 


4z/SMAMSW + LSND 




Ki 


= 0.0001, 


\k 


2 | = 0.002 


(0.75, 0.54, 0.99) 


(-0.04, 0.31) 


5z/SMA MSW + LSND 




Ki 


= |«2 


2 

) I 


K 2 \ 


= 0.032 


(0.88, 0.51, 0.96) 


(-0.09, 0.31) 



Wolfenstein parameters (p, 77) in some cases have significant corrections de- 
pending on the neutrino solution (see table [7]). In particular, for larger 
values of Ki, k 2 we obtain significantly larger negative values of p; in par- 
ticular consider p = —0.24 for the 3 v LMA MSW solution. This may be 
severely constrained by B — B mixing data. However in order to determine 
whether this is consistent with present data we must first extend our nu- 
merical analysis to include this process. We will look at this in a future 
paper |24[ . 



3 Neutrino Masses and Mixing Angles 

The parameters in the Dirac Yukawa matrix for neutrinos (eqn. |j) mixing 
v — v are now fixed. Of course, neutrino masses are much too large and we 
need to invoke the GRSY |2£| see-saw mechanism. 

Since the 16 of SO(10) contains the "right-handed" neutrinos v, one 
possibility is to obtain v — v Majorana masses via higher dimension operators 
of the form Q 

16 16 3 16 16 3 , ^ 16 16 3 16 16« a , 16 16 16 16 6 S ab . (9) 

The second possibility, which we follow, is to introduce SO (10) singlet 
fields iV and obtain effective mass terms V — N and N — N using only di- 

12 This possibility has been considered in the paper by Carone and Hall p7[ . 
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mension four operators in the superspace potential. To do this, we add three 
new SO(10) singlets {N a , a = 1,2; N 3 } with U(l) charges { -1/2, +1/2 }. 
These then contribute to the superspace potential 




where the field 16 with U(l) charge —1/2 is assumed to get a vev in the 
"right-handed" neutrino direction. Note, this vev is also needed to break the 
rank of SO(10). 

Finally we allow for the possibility of adding a U(2) doublet of SO(10) 
singlets N a or a U(2) singlet N 3 . They enter the superspace potential as 
follows - 



The dimensionful parameters //, /i 3 are assumed to be of order the weak 
scale. The notation is suggestive of the similarity between these terms and the 
/j, term in the Higgs sector. In both cases, we are adding supersymmetric mass 
terms and in both cases, we need some mechanism to keep these dimensionful 
parameters small compared to the Planck scale. 
We define the 3x3 matrix 



The matrix jl determines the number of coupled sterile neutrinos, i.e. there 
are 4 cases labeled by the number of neutrinos (N v = 3, 4, 5, 6): 

• (N v = 3) 3 active (// = /i 3 = 0); 

• (N u = 4) 3 active + 1 sterile 



WD fi'N a N a + fi 3 N 3 N 3 



(11) 



li = 



n' \ 
fj,' 

/1 3 J 



(12) 



V 



V 
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In this paper we consider the cases N v = 3, 4 and 5. 

The generalized neutrino mass matrix is then given by 



/ 


1 

m 

V o 



N v 

m 



ft V T 



N) 

\ 

V 
M N J 



(13) 



where 



and 



m 



Yp -j= sin/3 



(14) 



V 



3 K X e V 16 
(e 1 - 3 k 2 e) V w 




M 



N 



e' + 3 k 2 e) Vi e 
3eV 16 
re(l-a)T p ^ 16 

/ Ki S K 2 S 

k 2 S S 4> 
06 






K 6 



(15) 



Vi6; ^i6 are proportional to the vev of 16 (with different implicit Yukawa 
couplings) and S, <fi are up to couplings the vevs of S 22 , 2 , respectively. 

Since both V and are of order the GUT scale, the states v, N may 
be integrated out of the effective low energy theory. In this case, the effective 
neutrino mass matrix is given (at Mq) by Q (the matrix is written in the 



N) flavor basis where charged lepton masses are diagonal) 



U 



m (V 1 



M N V 



m 



-ji T v- 



m 



-m (v T y 





(16) 



13 This is similar to the double see-saw mechanism suggested by Mohapatra and 
Valle @. 

14 In fact, at the GUT scale Mq we define an effective dimension 5 supersymmetric 
neutrino mass operator where the Higgs vev is replaced by the Higgs doublet H„ coupled 
to the entire lepton doublet. This effective operator is then renormalized using one-loop 
renormalization group equations to Mz- It is only then that H u is replaced by its vev. 
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with 



e 



U e -- 



"a 



U e \ 
l) 



(17) 



U e is the 3x3 unitary matrix for left-handed leptons needed to diagonalize Y e 
(eqn. [|) and e , vq (e, v) represent the three families of left-handed leptons 
in the charged-weak ( -mass) eigenstate basis. 

The neutrino mass matrix is diagonalized by a unitary matrix U = U a f, 



rn 



diag = JJT 



(18) 



where a = {u e , u^, u T , u Sl , u S2 , u S3 } is the flavor index and i — {1, • • • , 6} is 
the neutrino mass eigenstate index. U a i is observable in neutrino oscillation 
experiments. In particular, the probability for the flavor state v a with energy 
E to oscillate into z/g after travelling a distance L is given by 



P{v a -> vp) = 5 af3 - 4 ^ U ak Up k U*j Upj sin 2 A jk 

k<j 



(19) 



where Aji 



and Sm % 



For N v < 4 we have 



1 1 
m j ~ m k- 



m,, 



m! Ill 



K\UJ ( K2 00 ( 



K\ e e r a 



b 

C23 
— u c 



C23 
C33 

-fc 



\ 

— u c 
-fc 

/ 



(20) 



where 
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C= (fft) (21) 

e 2 = (e') 2 + 9 («i - k|) e 2 

r< _ ^1 i 3 Kl r 2 a (3-4er) \ 3 (ki -k 2 .) e 2 -K 2 e e' j- 

O23- + 2 P ) P raC, 

Q _ 6ki <? to ( _i_ 3 Ki e 3 r 2 cr (1-cr) \ , ki e 2 r 2 a 2 /■ 



M = r cr e 

M3 Vl6 

a; mi 



C 



and 



y _|_ 3e 3 Kir 2 o-(l— cr) 



, A 2 v 2 sin 2 8 uid> m+ oj d> . . 

^16 ^16 ^16 ^16 

where in the approximation for mf we use 

f 

mt(= m top ) ~ A — ^ sin/3, (23) 

valid at the weak scale. 

In addition, for N u = 5 the off-diagonal piece of the mass matrix in eq . (|i~G|) 

reads 

-*™=-^(; ; in^O- (24) 



with 



d = ^ (25) 
m* (p 

g = (3eK 2 + e')/f (26) 

h = -3eKi/e 2 . (27) 



3.1 Three Neutrinos 

Before we discuss the case with non-zero K(i,2), let's recall the problem when 
K (i,2) — 0. For three active neutrinos with minimal family breaking vevs, 
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(0 2 ), (S 22 ), (A 12 ) ^ and k x = k 2 = 0, we find (at M G ) in the (u e , v„, u T ) 
basis 



to,. 



m' UT 



( 

6 1 I E/e (28) 



\ o i o 

is given in terms of two independent parameters {to', 6 } (see equa- 



tions [H], P2"|) . Note, this theory in principle solves two problems associated 
with neutrino masses. It naturally has small mixing between v e — since 
the mixing angle comes purely from diagonalizing the charged lepton mass 
matrix which, like quarks, has small mixing angles. While, for b < 1, — v T 
mixing is large without fine tuning. Also note, in this theory one neutrino 
(predominantly u e ) is massless. 

We have checked that in this theory it is possible to simultaneously fit 
both atmospheric and LSND data. We however cannot simultaneously fit 
both solar and atmospheric neutrino data. As discussed in paper I [14| this 



problem can be solved at the expense of adding a new family symmetry 
breaking vev 

(0 1 ) = M0 2 )- (29) 

In this paper we consider the most general family symmetry breaking 
vevs, given in equation [[], introducing two new complex parameters K\, k 2 . 
This allows us to obtain a small mass difference between the first and second 
mass eigenvalues which was unattainable before in the large mixing limit for 
v n — v T . Hence good fits to both solar and atmospheric neutrino data can 
now be found. In addition, in the previous section we showed that small 
values of «i 2 are consistent with good fits for charged fermion masses and 
mixing angles. In the next section we discuss these new solutions. 

4 Neutrino oscillations [ 3 active only ] 

In this section we consider the solutions to atmospheric and solar neutrino os- 



cillations with three neutrinos. The mass matrix is given in equation 20 with 



the parameter c = 0. There are three possible solutions to the solar neutrino 



15 This additional vev was necessary in the analysis of Carone and Hall. 27 
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data defined as small mixing angle [SMA] MSW, large mixing angle [LMA] 
MSW or "Just so" vacuum oscillations ||. In all three cases atmospheric 
neutrino data are predominantly described by — > v T oscillations. 

Instead of fitting the data directly, we compare our models with existing 2 
neutrino oscillation fits to the data |J . We use the latest two neutrino fits to 
the most recent Super-Kamiokande data for atmospheric neutrino oscillations 
and the best fits to solar neutrino data including the possibility of "just so" 
vacuum oscillations or both large and small angle MSW oscillations 0, [I], §]. 

For atmospheric neutrino oscillations we have evaluated the probabil- 
ities (P(Vfi — > Vp), P(v^ — > u x ) with x = {e, r, s}) as a function of 
x = Log[(L/km)/(E/GeV)]. In order to smooth out the oscillations we have 
averaged the result over a bin size, Ax = 0.5. In figures 2a and 4a we see that 
our results are in good agreement with the values of 5m 2 atm and sin 2 29 atm as 
given. 

For solar neutrinos we plot, in figures 3(a,b) and 5(a,b), the probabilities 
(P{y e v e ), P{y e — » u x ) with x = {/x, r, s}) for neutrinos produced at the 
center of the sun to propagate to the surface (and then without change to 
earth), as a function of the neutrino energy E u (MeV). We then compare 
our model to a 2 neutrino oscillation model with the given parameters. 



4.1 3 v SMA MSW solution 

In tables |8| and |9| we give the parameters for the fit corresponding to figures 
2(a,b) and 3(a,b). This model is indistinguishable from the results of the 
given parameters for 2 neutrino oscillations z/ M — v T for atmospheric data and 
v e — ^active for solar data. 

In order to obtain a SMA MSW solution we need to choose K\ = k\ 
to high accuracy. Note this value of «(i,2) corresponds to the only solution 
obtained previously (in I) with non zero k defined by ((f) 1 ) = n(4> 2 )- In 
fact, an SU(2) rotation of this case to zero gives non zero (S" 11 ), (S* 12 ) 
satisfying the relation k\ = k\. 

The parameter m' is determined by the high see-saw scale. Given m'(eqn. 
22] and table we find Vi§V{ & /<f> = 1.33 • 10 16 GeV which is consistent with 



the GUT scale. The large value of b (eqn. |2Tp results from S ~ 10 which 

16 For this calculation use an analytic approximation necessary to account for both large 
and small oscillation scales. For the details, see the appendix. 
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Table 8: Fit to atmospheric and solar neutrino oscillations [3 v 
SMA MSW] 

Initial parameters: ( K\ — /cf, | ac 2 I = 0-028 ) 
m! = 3.35 -lO" 3 eV , b = 15.0, $ fc = 3.30rad 



Observable 


Computed value 


5r n 2 atm 
Sill 2 29 aim 

sin 2 29 soi 


3.5 ■ 10" 3 eV 2 

0.99 
6.3 • 10~ 6 eV 2 
5.2 • 10~ 3 



Table 9: Neutrino Masses and Mixings [3 v SMA MSW] 

Mass eigenvalues [eV]: 0.000001, 0.0025, 0.059 

Magnitude of neutrino mixing matrix U a j 
i — 1, • • - , 3 - labels mass eigenstates. 
a = {e, /i, t} labels flavor eigenstates. 

0.997 0.0360 0.0599 " 
0.0677 0.672 0.738 
0.0172 0.740 0.673 



is needed in order to have one large and two small eigenvalues. 
4.2 3 v LMA MSW solution 

In figures 4(a,b) and 5(a,b) we present the comparison to a two neutrino 
oscillation model for atmospheric and solar neutrino data (see also tables [TO 
and [II]). For atmospheric data the fit is good for values of 



Log((L/km)/(E/GeV)) < 4 (see figures 4(a,b)) where the oscillations are 
predominantly given by v^ — v T . For larger x > 4 the probability Piv^ — > v^) 
is significantly smaller (~ 30 %) in our model than in a two neutrino model. 
This is due to the onset of significant z/ M — > v e oscillations. These larger 
values of x may be accessible in atmospheric oscillations. The maximum 
distance L for neutrinos of order 13,000 km, for upward going neutrinos, and 
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1.0 
J 0.8 
0.6 

CM 

0.4 
0.2 



0.0 



Atmospheric neutrino analysis (lie,! = 0.028, k 1 = k 2 ) 



♦ ♦ ♦ 



♦ ♦ ♦ ♦ 



♦ 3 neutrino model 
2 neutrino model with (sin ; 2e,8m 2 /eV ; ) = (0.99, 3.5x10 3 ) 



0.0 1.0 2.0 3.0 4.0 5.0 6.0 
Log[(L/km)/(E/GeV)] 



Figure 2 a: Probability Pivp, — > v^) for atmospheric neutrinos [3 v SMA 
MSW]. For this analysis, we neglect the matter effects. 



Atmospheric neutrino analysis (Ik,I = 0.028, K 2 = ) 



1.0 
> 0.8 
0.6 
0.4 
0.2 



0.0 



Ox = e 

♦ X = T 



♦ ♦ ♦ ♦ 



0.0 1.0 2.0 3.0 4.0 
Log[(L/km)/(E/GeV)] 



5.0 6.0 



Figure 2 b: Probabilities Piy^ 
neutrinos [3 v SMA MSW]. 



v x ) [x 



e, r and s) for atmospheric 
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Solar neutrino analysis (Ik,I = 0.028, K t = K 2 ") 



3 neutrino model 

2 neutrino model with (sin 2 28,8m 2 /eV 2 ) = (5.2x10 6.3x10 6 ) 




E v (MeV) 

Figure 3 a: Probability P[y e — > u e ) for solar neutrinos [3 v SMA MSW]. 



Solar neutrino analysis (iKjl = 0.028, K, = K 2 2 ) 




E v (MeV) 

Figure 3 b: Probabilities P[y e — > v x ) (x = fi, t and s) for solar neutrinos 
[3 v SMA MSW]. 
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Table 10: Fit to atmospheric and solar neutrino oscillations [3 v 
LMA MSW] 

Initial parameters: ( \k\ \ = 0.055, |k 2 | = 0.31 ) 
m! = 4.93 -10~ 2 eV , b = 0.84, $ 6 = -0.41 rad 



Observable 


Computed value 


5r n 2 atm 
Sill 2 29 aim 

sin 2 29 soi 


3.7 ■ 10~ 3 eV 2 

0.99 
2.3 • 10~ 5 eV 2 

0.77 



Table 11: Neutrino Masses and Mixings [3 v LMA MSW] 

Mass eigenvalues [eV]: 0.002, 0.005, 0.061 

Magnitude of neutrino mixing matrix U a j 
i — 1, • • - , 3 - labels mass eigenstates. 
a = {e, /i, r} labels flavor eigenstates. 

" 0.857 0.513 0.049 " 
0.368 0.563 0.740 
0.362 0.648 0.671 



the minimum detectable energy of order 0.1 GeV, corresponds to a value of 
xmax ~ 5. On the other hand, it would require a much more detailed analysis 
to determine whether our model is consistent with the data for fully contained 
events in the sub GeV ( < 1.33 GeV ) regime. We also note that this effect 



has been considered, in a recent analysis by Peres and Smirnov p9| , as a 
possible tool to distinguish the LMA MSW solution from the other solutions 
to the solar neutrino problem. 

A large mixing angle oscillation solution is obtained by tuning the lightest 
two neutrinos to be approximately degenerate with a near bi-maximal mixing 
matrix (see tables |H] and |TT| ), where the bi-maximal mixing matrix is given 
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Atmospheric neutrino analysis (hcj = 0.055, IkJ = 0.31) 



1.0 



>* 0.6 



♦ ♦ ♦ 

♦ 

♦ ♦ 

♦ 

♦ 




♦ 3 neutrino model 

y 2 neutrino model with (sin 2 28,Sm 2 /eV 2 ) = (0.99, 3.7xl0~ 3 ) 







0.0 1.0 2.0 3.0 4.0 5.0 6.0 
Log[(L/km)/(E/GeV)] 

Figure 4 a: Probability Pivp, — > v^) for atmospheric neutrinos [3 v LMA 
MSW]. For this analysis, we neglect the matter effects. 



Atmospheric neutrino analysis (IkJ = 0.055, Ik,I = 0.31) 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



o x = e 
♦ x = % 
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0.0 1.0 2.0 3.0 4.0 5.0 6.0 
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Figure 4 b: Probabilities Piy^ 
neutrinos [3 v LMA MSW]. 



v x ) [x 



e, r and s) for atmospheric 
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Solar neutrino analysis (Ik,I = 0.055, Ik,I = 0.31) 



| 0.6 
0.4 
0.2 
0.0 







♦ ♦ 3 neutrino model 

2 neutrino model with (sin 2 2e,8m 2 /eV 2 ) = (0.77, 2.3x10 5 ) 


^♦>MW * — * * * ♦♦«♦< 









E v (MeV) 

Figure 5 a: Probability P{y e — > u e ) for solar neutrinos [3 v LMA MSW]. 



1.0 

k 0.8 
f 0.6 

M 

0.4 
0.2 



0.0 l 
10 



Solar neutrino analysis (Ik,I = 0.055, Ik 2 I = 0.31) 



> x = n 

> X = X 



1 10 

E v (MeV) 



10' 



Figure 5 b: Probabilities P[y e — > v x ) (x = fi, t and s) for solar neutrinos 
[3 v LMA MSW] . 
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by H 



\U a 



0.71 
0.5 
0.5 



0.71 0.0 
0.5 0.71 
0.5 0.71 



(30) 



Note, a major difference in our case is the non-zero value for U Ve 3 ~ .049. 
However the constraint U Ve $ = chosen to satisfy CHOOZ M is much too 
strong. It is easy to see that our model is consistent with the null results of 
CHOOZ, i.e. 



V,- 



l-4|^ e3 | 2 (l-|^ 



3 | 2 )sm 2 



■ 1.275m 2 atm (eV 2 )L(km) < 
E v (GeV) 



> 0.98 



(31) 



for values of \U Ue s\ < 0.16 JJTJ. Finally the parameter b ~ 1 requires no fine 
tuning and given m' we find the high see-saw scale Vi6VY 6 /</> = 8.78 ■ 10 14 
GeV. 



4.3 3 v "Just So" Vacuum solution 

A vacuum solution is obtained by tuning the lightest two neutrinos to be 
even more degenerate than in the previous LMA MS W case. Our results are 



given in tables [12| and [13|. We have not given any figures since the results 



are standard vacuum oscillations. Once again we obtain a near bi-maximal 



mixing matrix |30[ with however U Ve % ~ .049. Nevertheless this model is 
consistent with CHOOZ data [Q (see the discussion of this in the LMA 
MSW case). Finally given the overall scale m' we determine the high energy 
scale to be V 1& VlJ<j) = 2.27 • 10 15 GeV and 6-1. 

In the next section we discuss a four neutrino solution to atmospheric, 
solar and LSND neutrino data in the theory with «(i,2) 7^ 0. 



5 Neutrino oscillations [3 active + 1 sterile] 

In the four neutrino case the mass matrix (at Mq) is given by equation |20| 
with c 7^ 0. 

As in the previous case of three neutrinos, we compare our model with 
two-neutrino oscillation models which have already been fit to the data [|1], 
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Table 12: Fit to atmospheric and solar neutrino oscillations [3 v 
Vacuum] 

Initial parameters: ( |/ci| = 0.004, | ac 2 I = 0.025 ) 
rri = 2.92 -lO" 2 eV , b = 1.73, $ 6 = -0.33 rad 



Observable 


Computed value 


^ m atm 

sin 2 29 atm 

5m 2 sol 
sin 2 26 sol 


3.5 • 10~ 3 eV 2 

0.99 
7.9 ■ 10- 11 eV 2 

0.97 



Table 13: Neutrino Masses and Mixings [3 v Vacuum] 

Mass eigenvalues [eV]: 0.00106037, 0.00106041,0.059 

Magnitude of neutrino mixing matrix U a j 
% — 1, • • • , 3 - labels mass eigenstates. 
a = {e, fi, t} labels flavor eigenstates. 

" 0.759 0.649 0.049 " 
0.429 0.513 0.744 
0.489 0.563 0.667 
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Table 14: Fit to atmospheric, solar and LSND neutrino oscillations 
[4 neutrinos SMA MSW + LSND] 

Initial parameters: = 0.0001, \kz\ — 0.002 

m! = 0.979 eV , b = -0.054, c = 0.101, $ fe = 5.59rad 



Observable 


Computed value 


^ m atm 


3.5 ■ 10~ 3 eV 2 


Sin 2 29 aim 


1.0 




5.0 ■ 10~ 6 eV 2 


sin 2 26 so i 


3.0 • 10~ 3 


^ m LSND 


0.53 


sin 2 29 LSND 


0.018 



|^, |5| . The results for our best fit are found in tables 
6(a,b), 7(a,b) and 8. 

In fig. 6a we evaluate P{y^ — > v^) where we also include a multiplica- 
tive fudge factor a. This is justified by the theoretical uncertainty in the 
normalization of the incident flux. Recall the observed number of muon 
neutrinos is given by 

N{v ll )= N (v„) P(iv-> ^) (32) 

where Nq{v^) is the theoretically expected incident neutrino flux which has 
an uncertainty of order 20%. We let N^v^) = Nf it ■ a where Nf it is the value 
used for the neutrino flux when fitting the data. 

We see that our result is in good agreement with the values of 5m 2 atm = 
3.5 ■ 10~ 3 eV 2 and sin 2 20 atm = 1.0 with a = 1.04. 

In fig. 6b we see that the atmospheric neutrino deficit is predominantly 
due to the maximal mixing between — v T , as in the case with k\ = k,2 = 0. 
However, in the case with K\ — k-i = there is also a significant (~ 10% 
effect) oscillation of — v s . In this case, that effect has vanished. This also 
means that sterile neutrinos will not come into thermal equilibrium in the 
early universe, due to the small mixing angle. Hence, at the nucleosynthesis 
epoch this model has only three neutrino species in thermal equilibrium. 

For solar neutrinos we see in fig. 7a that our model reproduces the neu- 
trino results for Sm 2 sol = 5m\ 2 = 5 x 10~ 6 eV 2 and a 2 neutrino mixing angle 



14 and hq and figures 
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Table 15: Neutrino Masses and Mixings [4 neutrinos SMA MSW 
+ LSND] 

Mass eigenvalues [eV]: 0.00002, 0.0022, 0.7248, 0.7272 
Magnitude of neutrino mixing matrix U ai 

i — 1, • • - , 4 - labels mass eigenstates. 

a = {e, /i, r, s} labels flavor eigenstates. 

0.997 0.0254 0.0480 0.0482 " 

0.0703 0.1079 0.7022 0.7003 

0.273 • 10~ 3 0.0292 0.7053 0.7083 

0.0181 0.9934 0.0852 0.0745 



Atmospheric neutrino analysis (Ik,I = 0.0001, k 2 l = 0.002) 



1.0 

- 0.8 

> 

L 0.6 

> 

s 0.4 
0.2 



0.0 
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♦ 4 neutrino model (a = 1.04) 
2 neutrino model with (sin 2 29.8m 2 /eV 2 , a) = (1, 3.5x10 3 , 1) 



0.0 1.0 2.0 3.0 4.0 5.0 6.0 
Log[(L/km)/(E/GeV)] 



Figure 6 a: Probability Pivp, — > v^) for atmospheric neutrinos multiplied 
by a, a fudge factor introduced to account for the uncertainty in the normal- 
ization of the incident flux. For this analysis, we neglect matter effects. 
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Atmospheric neutrino analysis (Ik,I = 0.0001, k 2 l = 0.002) 
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6.0 



Figure 6 b: Probabilities Piv^ 
neutrinos 



u x ) (x = e, r and s) for atmospheric 



sm 2 29 so i = 3 x 10 3 . The solar neutrino deficit is predominantly due to the 
small mixing angle MSW solution for v e — v s oscillations. The results are 



summarized in tables 14 and 15 



A naive definition of the effective solar mixing angle is given by 

sin 2 2fl 12 = 4 \\U el \\ 2 \\U e2 \\ 2 . (33) 

We note that the naive definition of sin 2 2#i2 underestimates the value of the 
effective 2 neutrino mixing angle. The fit value corresponds to sin 2 26 , 12 = 
2.6 x 1(T 3 . 

In fig 7b we see that oscillations into any active neutrino is substantially 
suppressed. This is unlike the case with K\ = k 2 = where there is also a 
significant (~ 8%) probability for v e — > v^. 

Finally with non vanishing k 12 we are now able to simultaneously fit 
atmospheric, solar and LSND data. This result is shown in figure 8 where we 
plot the probability Piv^ — ► v e ) as a function of neutrino energy relevant for 
LSND for our model compared to a two neutrino model with sin 2 29 = 0.018 
and 5m 2 = 0.53 eV 2 in the LSND allowed region Q This is in contrast 

17 Note, the probability for — > v e oscillations is almost identical. 
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Solar neutrino analysis (IkJ = 0.0001, Ik 2 I = 0.002) 



4 neutrino model 



2 neutrino model with (sin 2 20,8m 2 /eV 2 ) = (3x10 , 5xl0" 6 ) 




E v (MeV) 

Figure 7 a: Probability P[y e — ► v e ) for solar neutrinos 



Solar neutrino analysis (IkJ = 0.0001, Ik 2 I = 0.002) 
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0.006 



0.004 



0.002 



LSND neutrino analysis (IkJ = 0.0001, Ik,I = 0.002) , L = 30 m 




4 neutrino model 



0.000 



30 35 40 45 50 55 60 
E v (MeV) 



Figure 8: Probability Piy^ — > v e ) for LSND energies. 



to the case «i 2 = (paper I) where this was not possible. 

We now consider whether the parameters necessary for the fit make sense. 
We have three arbitrary parameters. We have taken b and c complex, while 
any phase for m' is unobservable. A large mixing angle for — v T oscillations 

is obtained with |6| ~ 0.05 [table [14[]. This does not require any fine tuning; it 

s v f 

is consistent with , ^jj ~ 0.17 which, taking into account Yukawa couplings, 
is perfectly natural (see eqn. |21|). The parameter c [eqn. and table |14j 
m 0.10 ~ H^m^tf) i m phes ^3 ~ 26(-^-) GeV. Considering that the standard 
/i parameter (see the parameter list in the captions to table ^ with value 
[i = 110 GeV and /x 3 [eqn. |ll| may have similar origins, both generated 
once SUSY is spontaneously broken, we feel that it is natural to have a light 
sterile neutrino. Lastly consider the overall scale of symmetry breaking, i.e. 
the see-saw scale. We have m' = 0.979 eV [table [14]] ~ ^ w y f . Thus we 

find Vw , V * 6 ~ m * f ~ 4.66 x 10 13 GeV. This is admittedly somewhat small 
but perfectly reasonable for (16) ~ (<f) 2 ) ~ Mq once the implicit Yukawa 
couplings are taken into account. 
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Table 16: Fit to atmospheric, solar and LSND neutrino oscillations 
[5 neutrinos SMA MSW + LSND] 

Initial parameters: = | ac 2 1 2 ; \^2\ — 0.032 

ml = 0.8380 eV , b = 0.9015, d = 0.0016, $ 6 = -3.18rad, $ d = -4.83rad 



Observable 


Computed value 


^ m atm 


3.7- 10~ 3 eV 2 


Sin 2 29 aim 


0.99 




5.7- 10~ 6 eV 2 


sin 2 26 so i 


4.0 • 10~ 3 


^ m LSND 


0.36 


sin 2 29 LSND 


0.026 



6 Neutrino oscillations [3 active + 2 sterile] 



In this case we have // ^ 0, /x 3 = (see eqn. [12]) where // sets the scale for 
the terms (mj^, (m„)5 a for a = {e, /x, r} (eqn. ^4|). We are able to find 
a good solution to atmospheric neutrino oscillations with maximal — > i/ r 
mixing, a solution to solar neutrino oscillations in the SMA MSW region and 



a fit to LSND. The fit is presented in tables [16] and [17| and in figures 9(a,b), 
10(a,b) and 11. 

Note, the parameter d (table [16| and eqn. ^4|) = Thus //' = 

m t d^r = 0.28 t^- GeV. In addition, we have m' = 0.838 eV [table M 

V 16 V 16 ' 1 

« Thus we find V i6 ~ ~ 5.3 x 10 13 4- GeV. In order 

to obtain this solution without fine tuning we must assume that the ratio 
T/T~ ~ 100. As in the previous four neutrino case, this may be attributable 

M.6 

to ratios of Yukawa couplings. 



7 Discussion 

In this paper we analyze the predictions for both charged fermion and neu- 
trino masses and mixing angles in an SO(10) SUSY GUT with U(2)xU(l) n 
family symmetry. We find that, if we allow for the most general family 
symmetry breaking vevs, the model can accommodate three different three- 
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Table 17: Neutrino Masses and Mixings [5 neutrinos SMA MSW 
+ LSND] 



Mass eigenvalues [eV]: 0.88 • 10~ 7 , 0.0007, 0.0025, 
Magnitude of neutrino mixing matrix U a ; 

i = 1, • • • , 5 - labels mass eigenstates. 

a = {e, /i, r, Si, s 2 } labels flavor eigenstates. 



0.6013, 0.6043 



0.0586 0.9940 0.0297 

0.0033 0.0802 0.0182 

0.0018 0.0356 0.0617 

0.0036 0.0291 0.9975 

0.9983 0.0585 0.0053 



0.0763 0.0430 

0.6998 0.7096 

0.7091 0.7015 

0.0401 0.0507 

0.0014 0.0015 



Atmospheric neutrino analysis (IkJ = 0.001024, Ik,I = 0.032) 
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Figure 9 a: Probability Pivp, — > 
analysis, we neglect matter effects. 



for atmospheric neutrinos. For this 
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Atmospheric neutrino analysis (IkJ = 0.001024, Itql = 0.032) 
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Figure 9 b: Probabilities Piv^ — > z/r) (a? = e, r, s x and s^) for atmospheric 
neutrinos. 
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Solar neutrino analysis (Ik,I = 0.001024, = 0.032) 



x = \i 




E v (MeV) 

Figure 10 b: Probabilities P{y e — > v x ) (x = //, r, si and S2) for solar 
neutrinos. 



LSND neutrino analysis (Ik,I = 0.001024, licj = 0.032), L = 30m 

0.006 i 




5 neutrino model 

x 2 neutrino model with (sin ! 29.8m z /eV 2 ) = (0.0259, 0.363) 



o.ooo 

30 35 40 45 50 55 60 
E v (MeV) 

Figure 11: Probability Piy^ — > v e ) for LSND energies. 
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neutrino oscillation solutions to atmospheric and solar neutrino data, one 
four and one five neutrino solution to atmospheric, solar and LSND data. 
We also find a three neutrino solution to atmospheric and LSND data alone. 
In spite of all this freedom in the neutrino sector, the fits for charged fermion 
masses and mixing angles are relatively unaffected. 

In all cases we find atmospheric neutrino data described by maximal 
v n ~> v t mixing.^ Super-Kamiokande is able to distinguish — > v x for 
v x = v T or v s (see talks by K. Scholberg and W.A. Mann [fj]). There are two 
proposed methods. The first uses the measured zenith angle dependence, 
since there is an MSW effect in the earth for v s but not for v T . This effect 
suppresses — > v s oscillations for high energy neutrinos coming from below. 
Recent data does not show such an effect; thus favoring v x = v T . The 
second method uses the ratio of neutral current [NC] to charged current [CC] 
processes which can distinguish between the two. Here there is preliminary 
data favoring v x = v T . This ratio satisfies 

R(nc/cc) < 1 for z/ x = v s (34) 
= 1 for u x = v T . 

Using SuperK data for n° events produced by neutral current neutrino scat- 
tering in the detector one measures 



E> — Qi~"/e) CQiQ /oe\ 

K (NC/CC) = W/e) MmteCarlo ^ 

= 1.11 ± 0.06(data stat.) ± 0.02(MC stat.) ± 0.26(sys.) 

The oscillations —>■ v T may also be visible at long baseline neutrino 
experiments. Both K2K and MINOS JlTJ are designed to test for 



disappearance. For example at K2K ||10|| , the mean neutrino energy E 



1.4GeV and distance L = 250 km corresponds to a value of x = 2.3 (see 
figures 2a, 4a, 6a and 9a) and hence Piy^ — > v^) ~ .45. 

Results on solar neutrino oscillations or LSND will, on the other hand, be 
able to narrow down the acceptable regions of parameter space, but cannot 
test this class of models. 



18 We have not searched for solutions with maximal — * v s mixing, since this is not 
favored by the latest Super-Kamiokande data M. 
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Finally, the plethora of solutions presented in this paper is in stark con- 
trast to the unique solution obtained assuming the minimal family symmetry 
breaking vevs studied previously in paper I |14[| . In the latter case we cannot 
find any three family solutions to both atmospheric and solar data and we 
find a unique four neutrino solution to atmospheric and solar data but NOT 
LSND. Thus it is clear that the neutrino sector is in general much less con- 
strained than charged fermions. Nevertheless, it is pleasing to find a simple 
SUSY GUT which can accommodate all of this low energy data. 
Acknowledgements This work is partially supported by DOE grant DOE/ER/01545- 
767. 
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8 Appendix 



Solar neutrino analysis 

In this appendix we describe in detail the approximation which we used 
in the numerical analysis of solar neutrino oscillations. The Schrodinger 
equation for solar neutrinos is given by 

4*?(*) = H a ^ u (t), (36) 



dt 



H aP = ^r!:^ + VaWap. (37) 

Ah/ 



Here ty" is a state vector for neutrinos with flavor a (a = e, fi, r, and s for 
four neutrino model 0), H is the Hamiltonian for solar neutrinos, and E is 
the neutrino energy. The mass matrix m u in the flavor basis is given by (see 



equation [18]) 



m 



v 

flavor 



where U is the mixing matrix for neutrinos (v^ vor = U Q 
for four neutrino model) and m d J a9 is the diagonal mass matrix in the mass 
eigenstate basis. V a {t) is a time-dependent potential for neutrinos with flavor 
a as follows: 

V e (t) = V2G F {n e (t)-~n n (t)}, 

V,(t) = V T (t) = -V2G F ^n n (t), 

V,(t) = 0, " (39) 

where Gp is the Fermi coupling constant. Here we assume that electron (n e ) 
and neutron (n n ) number densities at a distance r = ct from the center of 
the sun are given by 

n e = 4.6 x 10 n exp (^-10.5^ eV 3 , (40) 

n n = 2.2 x 10 n exp (^-10.5^ eV 3 , (41) 



19 Here we present our method of solar neutrino analysis in a four neutrino model. The 
method can be easily extended to a three neutrino model or a model with more neutrinos. 
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where R is a solar radius. 

Mass scales for the atmospheric and LSND neutrino problems {Sm 2 tm ~ 
1CT 3 eV 2 , <5m 2 SND ~ 1 eV 2 ) are much larger that for the solar neutrino prob- 
lem (5m 2 olar < 1CT 5 eV 2 ). When we include the mass scales for atmospheric 
and/or LSND neutrinos and solve the Schrodinger equation for the solar 
neutrino problem, it is almost impossible to solve it numerically because of 
these larger mass scales and the rapid fluctuations they produce. Thus, in 
order to solve the Schrodinger equation numerically, we use the following 
approximation. 

We divide the mass term m\m v into two parts: 



= Urn 2 U ] + U m 2 



2 Tj\ 
H u i 



(42) 



where m L (m 2 H ) is a "Light" ("Heavy") part, 



m T = 



m 



H 



( m\ 



( o 



V 



/ 



m 2 ) 



(43) 



(44) 



and we assume that 5m 2 x 



m\ is the scale for solar neutrino problem 



and m\ < m 2 , <^ < m\. Then the Hamiltonian H is given as 



H = H H + H L , 



H l = U^U^ + V a (t)5 aP . 
The state vector is also divided into two parts as follows: 



(45) 
(46) 

(47) 
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where we define A to satisfy the following equation: 

. d 



H H A(t), 
I 



A(t = 0) 

where / is a unit matrix. We can easily solve the equation 
solution is given by 

A(t) = exp(-iH H t), 



U exp 



m 



2E 



(48) 
and the 



(49) 



Then <&,, satisfies 



2E 



exp %■ 



2E J 



U exp 



.m 



H 



2E 



t 



U^ v (t). (50) 



Since the mass scales 7713,4 included in the matrix m 2 H are too large for 
MSW effects, the exponential terms exp (^hi^^-t^j oscillate rapidly There- 
fore we replace them by their time-averaged values: 



m 



exp ±,— t 







(51) 



Then equation [5(] has the following approximate form 



m 



Li 



2E 



y 7 i+2^i+2 



(52) 



where the indices a, (3, 7 run from 1 to 4, on the other hand, the indices 
i, j from 1 to 2. We then solve equation |52] with the initial condition 

$ u {t = 0) = (1,0,0,0) or #„(f = 0) = (1,0,0,0). (53) 

Finally, the oscillation probability P(u e — > u a ) (a = e,fi,r or s) at time 
t is given by 



P(v P 



E u ai ul 

1=1,2,^=1-4 



+ E 

8=3,4 



E 

13=1-4 



.(54) 
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where the ~ in the last line (equation |54[) refers to the fact that the time 
average of exp(±i-^rt) was used. 
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